Nonclassical-light generation in a photonic band-gap nonlinear planar waveguide 
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Optical parametric process occurring in a photonic band-gap planar waveguide is studied from 
the point of view of nonclassical-light generation. Nonlinearly interacting optical fields are described 
by the generalized superposition of coherent signals and noise using the method of operator linear 
corrections to a classical strong solution. Scattered backward-propagating fields are taken into 
account. Squeezed light as well as light with sub-Poissonian statistics can be obtained in two-mode 
fields under the specified conditions. 

PACS numbers: 42.50.-p Quantum optics, 42.65.-k Nonlinear optics, 42.50.Dv Nonclassical states of the 
electromagnetic field, 42.65.Yj Optical parametric oscillators and amplifiers 



I. INTRODUCTION 

During the last years an increasing attention has been 
devoted to properties of photonic band-gap structures 
PIQ- It has been shown that intense nonlinear processes 
can occur in such structures built up from nonlinear ma- 
terials. For example, second harmonic and sub-harmonic 
generation has been predicted and also observed 0, • 
Advantages of these structures from the point of view of 
nonlinear interactions are based on high densities of local 
optical modes, on spatial localization of optical modes in 
confined regions of the structure and in convenient ful- 
filling of phase-matching conditions of a given nonlinear 
process. Although quantum nonlinear optics in photonic 
crystals is promising, investigation of these structures has 
been performed mostly in classical domain. Such struc- 
tures are usually very short (typically tens of /im) and 
so they can be conveniently used in microoptoelectron- 
ics as sources of light in near future. These nonlinear 
structures can also serve as sources of light with non- 
classical properties, as predicted in As discussed in 
, second-harmonic-generation process in a planar non- 
linear waveguide with a corrugation on the top of the 
waveguide can be used to control squeezing of the funda- 
mental field; the corrugation reproduces a photonic-band 
gap structure. Periodicity of the grating was selected to 
give rise to a longitudinal confinement of the pump field, 
phase matching of the process was achieved introducing 
a spatial modulation of nonlinear susceptibility. 

The aim of the paper is to investigate nonclassical 
properties of light (squeezing of vacuum fluctuations and 
sub-Poissonian photon-number statistics) generated in 
a nonlinear planar photonic-band-gap waveguide fabri- 
cated in such a way that a strong optical parametric pro- 



cess occurs inside. The corrugation of the waveguide is 
suitably selected in order to confine longitudinally signal 
and idler modes, while a weak longitudinal confinement 
of the pump field is assumed. Phase matching for the 
parametric process is achieved due to the presence of the 
grating and due to modal dispersion of the guided wave 
geometry (see, e.g., in 0), although our considerations 
are not specific to a particular material. 

The paper is organized as follows. Equations describ- 
ing nonlinear interaction of optical fields are derived in 
Sec. 2. Possibility of squeezed-light generation is dis- 
cussed in Sec. 3. Sec. 4 is devoted to generation of light 
with sub-Poissonian photon-number statistics. Sec. 5 
provides conclusions. Quantum derivation of the corre- 
sponding equations is contained in Appendix A, whereas 
Appendix B deals with commutation relations among op- 
erators of the outgoing fields. 



II. DERIVATION OF EQUATIONS 
DESCRIBING THE NONLINEAR INTERACTION 

An optical field present inside a planar waveguide 
can be decomposed into modes. Then we have for the 
electric-field amplitude E(r, t) at point r = (x. y, z) and 
time t the following expression: 
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x [A m (z)f m (x,y)exp[i(k. m ) z z - iu m t] - h.c] , (1) 

where A m is amplitude of the m-th mode, f m 
means transverse eigenfunction of the m-th mode 
(A x y f m (x,y) = 0), e m stands for polarization vector, 
w m denotes frequency and k m is wavevector of the m-th 
mode. Mean permittivity of the waveguide is denoted as 
e r , €q stands for vacuum permittivity, h is reduced Planck 
constant and V volume of the structure. The symbol h.c. 
stands for hermitian-conjugated terms. 
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The electric-field amplitude E(r, t) fulfils the wave 
equation inside the waveguide: 



and 



V 2 E 



<9 2 E 



d 2 P nl 

dt 2 



(2) 



where [i denotes vacuum permeability and P n i describes 
nonlinear polarization of the medium. Relative permit- 
tivity e r (r) can be written as follows: 



e r (x, y, z) = e r (x, y) [1 + Ae r {x, y, z)} 



(3) 



Small variations of permeability described by Ae r (x, y, z) 
are related to the corrugation caused by a linear grating 
fabricated on the surface of the waveguide. It is useful 
to decompose permittivity variations Ae r as a function 
of spatial coordinate z into harmonic functions: 



Ae r (x,y, z) 



OO 
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0: 



(4) 



e q are coefficients of the decomposition. Polarization P n ] 
of the medium is determined using second-order suscep- 
tibility tensor \'- 



P„i = e oX ■ E • E; 



(5) 



the symbol ■ stands for tensorial multiplication. 

Substitution of the expression for E(r,t) in Eq. JTJ 
into the wave equation in Eq. @ and the assump- 



tion 



c)-A,„ 



< |(k m ) z ^^| (analog of slowly- varying 



envelope approximation in time domain to spatial evo- 
lution) results in the system of nonlinear differential 
equations for amplitudes A m of modes participating in 
the nonlinear interaction. We have six modes: signal 
forward-propagating mode (with amplitude A sp ), sig- 
nal backward-propagating mode (A SB ), idler forward- 
propagating mode (Ai P ), idler backward-propagating 
mode {Ai B ), pump forward-propagating mode (A PF ), 
and finally pump backward-propagating mode (A PB ) . 
The system of differential equations is written as follows: 



dA s 



dz 

dA lB 
dz 

dA SE 
dz 

dA iE 
dz 

dAp f 
dz 

dA PE 
dz 



iK s exp(-iS s z)A SB 

+ 2K F exp(iS F z)A PF A* F 

iKi exp(-iSiz)A iB 

+ 2K F exp(iS F z)A PF A* SF 

—iK* exp(iS s z)A SF 

- 2Kb exp(-iS B z)A PB A* E 

—iK* exp(iSiz)A iF 

- 2K B exp(-iS B z)A PB A* 1 

~2K F exp(-iS F z)A SF A lF 
2K* B exp(iS B z)A SB A iB 



(6) 



\(k aF ) z \ + \(k aB ) z \-S l 



s,i, 



c 27T 

5l = A? 
S h = l(k 



Pb)z 



-\(K b ) z \-\(KU b = F,B. (7) 



The linear coupling constants K s and Ki are given as: 
l(k„J ; 



K„ 



dxdy e t (x, y)f* (x, y)f aB (x, y), 



(8) 



The expressions for the nonlinear coupling constants K F 
and K B are: 



Ki, 



I fJ,hu p UJ s UJi 



J dxdy x ■ e p ■ e s 



x fp b ( x ,y)fs h ( x iy)fi h ( x >y)> 



F, B. (9) 



The same equations as in Eq. 10 can be derived using 
quantum description of optical fields. The quantum ver- 
sion of Eq. © uses operators of electric-field amplitudes 
instead of classical electric-field amplitudes (for details, 
see Appendix A). A general solution of the quantum vari- 
ant of the nonlinear system of operator equations in Eq. 
© can only be found using a completely numerical ap- 
proach. However, taking into account conditions of the 
nonlinear interaction in a real photonic-band-gap planar 
waveguide, we can apply the approximation of small lin- 
ear operator corrections to a classical strong solution for 
mean values. We can write in this case: 



A a — A a + 5A a , 
s F ,i F ,p F ,s B ,i B ,p B . 



(10) 



Mean amplitudes A a obey classical equations given in Eq. 
©, whereas the evolution of small operator corrections 
SA a is governed by the following equations: 



dSA s 

dz 
dSA t 

dz 
d5A Si 

dz 
dSA ti 

dz 
d8A. 



VF 



dz 
d5A n 



= K, S 5A S 
= ICiSAij 
= K*M a 
= K*5Ai 
= -/CI 



- K, F 
K, F 

- IC B 

- IC B 



A PF SA ip 



A PF 5A\ F 



A PB sa* b 



A Pb 5A\ e 



A* F 5A PF 



A* SF 5A PF 



A*5A r . 



KJA PB 



A SF SAi F 



A lF SA s 



dz 



= K\ 



A SB 5Ai 



A iB 5A s 



(11) 



The constants IC S , Ki, IC F , and KL B introduced in Eq. 
are defined as: 



JC a = iK a exp(~iS a z). 
KL F = 2K F exp(z(5p;z), 
K, B = 2K B exp(— iS B z). 



a = s, i, 



(12) 
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Assuming linear interaction being much stronger than 
the nonlinear one (KpA and Kg A are much lower than 
K s and Kf, A means an arbitrary classical amplitude), 
Eq. © can be solved analytically. We then have for the 
signal and idler modes: 



AaF{z) 



.S a z 



exp —i 



x [A aF (0) cos(A Q z) + C a sin(A Q z)] , 



A aB (z) = exp ( i— 



X, (t)i ( -^cos(A a2 ) + ^sin(A a z) 
+ C a ( -;^sin(A a z) - ^-cos(A Q z) 
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\Ka\\ 



a = s.i. 



(13) 



The symbol A SF (0) [Ai F (0)] stands for an incident signal- 
[idler-] field classical amplitude. Values of the constants 
C s and Ci are determined from the conditions A SB (L) = 
and A iB (L) — 0; L denotes the length of the waveguide. 

Classical amplitudes A PF and A PB of the pump modes 
inside the waveguide are given as: 



A Pb (z) = A Pb (0) 



i(A b 



± -■ — A-^y ( exp [-^ Ab - As 

va^a' + A^ (expH(A 6 - A s 
i{A b - A s + AjJ 

' 2K '^ e ' (expH(A 6 + A s 



± 



± 



i(A b + A s - At) 



± 



~ Ki:F '- :F > (exp[-i(A fe + A s 



i(A b + A S + Ai) 



Ai)z] 
Ai)z] 

Ai)2] 
A,;)zl 



F,B. 



-1) 
-1) 
-1) 

(14) 



The upper signs in Eq. (|14fl are appropriate for A pF 
whereas the lower signs stand for A PB . The constants 
occurring in Eq. (|14|) are given as: 



A F 
A B 

f 

<->a B 
f 

J a B 



S s Si 
= 6p+ 2 + X 
: 5 S Si 

= -Or - — - — 



~(A ap (0)-iC a ) 

^(A aF (0) + iC a ), 



2A a - 5 a 
4K a 

2A Q + Sg 

4K„ 



(A aF (0)-iC a ) 

(-A ap (0)-iC a ), a = s,i.(15) 



The constant A PF (0) characterizes the incident pump- 
field amplitude and value of the constant A pB (0) is de- 
termined from the condition A PB (L) = 0. 



Solution of the system of linear equations in Eq. Hll|) 
for operator corrections can be found numerically and 
expressed in the matrix form: 



SAF,out 

SAbm 



(Uff 
\Ubf 



U fb \ 
U BB J 



where 



SAfsu = 



5 At 



fSA SF (0)\ 
SX F (0) 

sA PF (o) 

\SAl(0)J 
fSX B (L)\ 
5A\ B {L) 
SX B (L) 
SAt(L) 
SA PB (L) 

\saUl)J 



SA 



_F,out 



8A 



B,out 



SA F ,in 
SAB.out 



f$A SF (L)\ 
SAt F (L) 
SA iF (L) 
SAl F (L) 
5A PF (L) 
\5Al F (L)J 
fSX B (0)\ 
SAt B (Q) 
SX B (0) 
SA\ B (0) 
SA PB (0) 

\sa; b (o)J 



(16) 



(17) 



The matrices Uff, Ufb, Ubf, and Ubb are determined 
by numerical solution of Eq. (|llfl . Precision of numerical 
solution of Eq. I|ll|) can be monitored using identities 
stemming from commutation relations among operators 
(see Appendix B). 

Input-output relations among linear operator correc- 
tions can be found solving Eq. (|16|) with respect to vec- 
tors SAp.out and S Ab, out- 



SAp i0ut 

SAs.out 



( Uff — UfbU bb Ubf UfbU b b 
V —U bb Ubf U 1 



BB 



SABAn 



(18) 



The output linear operator corrections contained in vec- 
tors 8AF,out and (>AB,out obey bosonic commutation re- 
lations provided that the input linear operator correc- 
tions in vectors <5AF,in and SAbah obey bosonic com- 
mutation relations. It has been shown in Q that this 
nontrivial property is fulfilled by any system described 
by a quadratic hamiltonian. 

Eq. 1|11|) can be solved also iteratively assuming 
weak linear (IC S L <C 1, )CiL <C 1) as well as nonlin- 
ear (1CfA Pf L <C 1, KbA Pb L <C 1) interactions among 
modes in the planar waveguide (weak-interaction approx- 
imation). The obtained expressions provide a useful in- 
formation about the behaviour of physical quantities of 
interest (see the next sections). 



III. SQUEEZED-LIGHT GENERATION 

Squeezing of fluctuations of an optical field below the 
level characterizing vacuum fluctuations can be indicated 
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by mean values of the variances of quadrature compo- 
nents qj [q~j = Aj + Aj, Aj stands for electric-field- 
amplitude operator of mode j] and pj [fij — — i(Aj — .At)] 
in mode j. Maximum amount of available squeezing mea- 
sured under a suitably chosen value of the local-oscillator 
phase in homodyne-measurement scheme is given by 
principal squeeze variance Aj |3 ■ 

The above-introduced quantities can be generalized to 
optical fields composed of two modes; we have qj k — 
qj + q k and pjk = pj + Pk for quadrature components of 
the field composed of modes j and k. An optical field 
described by a compound mode can be obtained in an 
output port of a beamsplitter that combines two single 
modes from its inputs. 

We assume that the interacting fields can be described 
in the framework of the generalized superposition of sig- 
nal and noise |nj (coherent states, squeezed states as 
well as noise can be considered) and then we have [TT| : 

((A<z,) 2 ) = l + 2[B j+ Re(Cj)], 
((Ap,) 2 ) = l + 2[B j -Re(C j )], 

Xi = l + 2[Bj-|C,-|], (19) 
{(Aq 3k ) 2 ) = 2 [1 + B s + B k - 2Re(D jk ) 
+Re(Cj +C k + 2D jk )} , 



{{Ap ]k ) 2 ) = 2 [1 + Bj 
-Re{Cj 
Xj k = 2 [l + Bj 



-B k -2Re(D jk ) 
C k + 2D jk )], 
-B k -2Re(D jk ) 



-\a+C k + 2D 



(20) 



Symbol ( ) denotes quantum statistical mean value. 
Quantities Bj, Cj, Dj k , and Dj k are defined as 



Bj = (AA'jAAj), 

Cj = ((AA,) 2 ), 

Dj k = (AAjAA k ), 

Dj k = -(AijAi fc ,) 



3 + h 



(21) 



AAj = Aj — (Aj). Quantities defined in Eq. (|21|l can 
be expressed in terms of the matrices u and v defined 
in Appendix B and using incident values of Bj^j± an d 
Cj,in,A related to antinormal ordering of field operators 
(for details, see [Tlf'): 

Bj,in,A = cosh 2 (rj) +n ch ,j, 



C 



— exp(z?9j) sinh(2rj). 



(22) 



Symbol rj denotes squeeze parameter of the incident j-th 
mode, dj means squeeze phase, and n c h.j stands for mean 
number of incident chaotic photons. Values of {(Aqj) 2 ), 
{(Apj) 2 ), and Aj less than one mean squeezing in single- 
mode case. Squeezed light is generated in a compound- 
mode (two-mode) case if values of ((Aqj k ) 2 ), {(Apj k ) 2 ), 
or Ajfe are less than two. 

When discussing properties of the nonlinearly inter- 
acting modes we use symmetry based on the exchange of 



signal and idler modes. Quantities and properties that 
can be derived from this symmetry are not mentioned 
explicitly. 

Discussion of squeezing is based on the investigation of 
principal squeeze variances because they give the maxi- 
mally allowed amount of squeezing reachable in an ex- 
periment. 



A. Weak-interaction approximation 

Assuming incident coherent states (or vacuum states) 
principal squeeze variances of single modes have the fol- 
lowing form in weak-interaction approximation: 



As F 1 ~h 2\Ip F | , 



X S a = 1 + 214 



X PF — 1, 

X PB = 1- 



(23) 



and 



dz tC F (z)A PF (z), 



dz JC B {z)A PB (z). 



(24) 



The expressions contained in Eq. (|23|l are correct up to 
the second power of constants K s , K i: K F , and K B - 
They show that principal squeeze variances of single 
modes are larger than one and thus no squeezing can 
occur in single modes. 

The following expressions for the compound modes are 
reached under the same assumptions: 



A., 



2 1+LL 



Xsis>.i-a — 2 [l + l/nci | 2 



\I PB \ 2 ], 
2\I PF \ 2 ] 



= 2 1-217, 



Pb I 



2|7 PB | 2 ] 



2 [1+17^-217^1 , 



Xs F ,p B — 2 [1 + \Ip F 

^S B ,PB = 2 [l + \Ip B 
X SB: PF = 2 [l + \lp B 



2|7„ 



A 



Pf,Pb 



= 2, 



(25) 



where 



L l,p F 



L PB,S 



dz / dz'K,*{z)JC F {z')A PF {z') 1 
Jo 

dz / dz'K. B (z)A PB (z)IC s (z'), 



dz / dz'IC F {z)A sp {z))C F {z')A PF (z'), 



dz / dz'lC B (z)A PB (z)]C B (z')A SB (z'). 
o Jo 



(26) 
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According to Eq. I|25|) squeezing can occur in compound 
modes (sp,ip) and (sb^b) owing to nonlinear process 
among the forward-propagating modes and that among 
the backward-propagating modes. Light in compound 
mode (sf,*s) can be squeezed if the linear coupling 
(describing scattering of light in the photonic-band-gap 
waveguide) is stronger than the nonlinear one. As is 
seen from the expression for \ SF ,i B in Eq. I|25[l squeez- 
ing originates in the nonlinear process and linear cou- 
pling between forward- and backward-propagating fields 
is inevitable for 'transfer of squeezing' into this mode. If 
classical amplitudes of the forward- (backward-) propa- 
gating signal mode are greater than those of the forward- 
backward-) propagating pump mode, squeezing can be 
reached in compound mode (sf,Pf) [(sb,Pb)]- 



B. Numerical analysis of squeezing 

Complete analysis of the behaviour of interacting 
modes can be reached only numerically. 

We assume a strong incident forward-propagating 
pump field and also nonzero incident signal and idler 
forward-propagating fields. Squeezed light cannot be 
generated in single modes in this case. However, com- 
pound modes (sf,*f)j (sb,!b), and (sf,ib) provide 
squeezed light at the output under suitably chosen values 
of the waveguide parameters. 

Values of parameters characterizing a real nonlinear 
planar photonic-band-gap waveguide are assumed to lie 
around the point given by L = 2 mm, K s = Ki — 
5 mm -1 , Kp — Kb = 5 x 10~ 7 mm -1 mV _1 , and 
A PF — 10 7 Vm" 1 . That is why we concentrate our atten- 
tion on quantitative analysis of behaviour of the waveg- 
uide in the vicinity of this point. We note that equality of 
values of linear coupling constants K s and Ki can barely 
be found in a real waveguide but real modes can be cho- 
sen in such a way that values of these constants are nearly 
the same. The quantities 1 mm and 10 6 Vm -1 are used 
as appropriate units in the following graphs. 

The role of length L of the waveguide to squeezed-light 
generation is revealed in Fig. ^ We would like to remind 
that any change of length L of the waveguide is accom- 
panied by a little change of the period of corrugation, in 
order to maintain the same mode profile and the same 
longitudinal confinement properties of the modes. Fix- 
ing these criteria, the modes (sF,ip) and (sb^b) are 
squeezed for any length L of the structure (see Figs. 
mi,c). However, the maximum value of squeezing is 
reached for values of length L around 0.5 and then there 
is saturation in values of squeezing. Squeezing of these 
modes around 20 % can be reached. The optimum value 
of length L lies around 0.3 for mode (sf,«s) composed 
of one forward- and one backward-propagating field, as 
is shown in Fig. Squeezing in mode (sf,Ib) can oc- 
cur owing to 'scattering of the already generated squeezed 
light' between forward- and backward-propagating fields. 
Saturation of the behaviour for larger values of length L 




L 

FIG. 1: Variances {(Aq) 2 ) (solid line with A) and ((Ap) 2 } 
(★) and principal squeeze variance A (❖) for mode (sf,if) 
(a), (sf,ib) (b), and (sb,ib) (c) in dependence on length 
L of the structure; Kf = Kb = 5 x 10~ 2 , K s = Ki = 5, 
5 S = 5, = 0, S F = 8 B = 0, A PF = 10, A SF = A lp = 0.1, 
A PB = A SB = Ai B = 0, incident coherent states are assumed. 
Backward-propagating fields exit the waveguide at z = 0. 



is typical for interactions containing both forward- and 
backward-propagating fields . 

Values of squeeze variances monotonously decrease 
with increasing values of nonlinear coefficients Kp and 
Kb and with increasing incident pump-pulse amplitudes 
A PF . This clearly shows that the origin of nonclassi- 
cal properties of outgoing fields lies in the nonlinear 
three-mode interaction. As the graph in Fig. [3 indi- 
cates, an increase of values of nonlinear coefficient K n i 
(K n i — Kp = Kb) and pump- pulse amplitude A PF two 
times with respect to the values characterizing the 'fixed 
working point' of the waveguide leads to values of squeez- 
ing around 50 % for mode (sf^f)- 

Linear coupling between forward- and backward- 
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FIG. 2: Principal squeeze variance A of mode (sf,if) in de- 
pendence on nonlinear coupling constant K n i (K n i = Kf = 
Kb and input pump- pulse amplitude A PF (0); L = 2 and val- 
ues of the other parameters are the same as in Fig. 1. 




FIG. 3: Principal squeeze variance A of mode (sf,if) in de- 
pendence on linear coupling constants K s and Ki\ L = 2 and 
values of the other parameters are the same as in Fig. 1. 



propagating fields (described by coupling constants K s 
and Ki) originates in scattering of forward-propagating 
fields in the photonic-band-gap waveguide and does not 
lead to any nonclassical behaviour. Moreover, a stronger 
linear coupling between signal modes and idler modes 
suppresses squeezing, as is demonstrated in Fig. It 
holds for mode (sf,if) that the larger the values of lin- 
ear coupling constants the larger the values of quantities 
characterizing squeezing (see Fig. [SJ- However, nonzero 
values of linear coupling constants are necessary for ob- 
taining squeezing in modes (spyig) and (se^b)- The 
reason is that linear coupling enables energy transfer to 
backward-propagating fields (being in vacuum states at 
the input) from the forward-propagating fields. 

Phase matching of all optical-field interactions occur- 
ring in the structure cannot be usually reached in a real 
structure. A typical influence of nonzero linear phase 
mismatch is demonstrated in Fig. ^describing behaviour 
of modes (sF,ip), (sf^b), and (ss,is). A nonzero 
value of linear signal-fields phase mismatch 6 S causes os- 



2.2 n 
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FIG. 4: Principal squeeze variance A for modes (sf,if) (solid 
line with *), (sf,ib) (A), and (sb,«b) (o) in dependence 
on linear phase mismatch S s ; L = 2 and values of the other 
parameters are the same as in Fig. 1. 



dilations in quantities characterizing squeezing. It effec- 
tively suppresses the influence of linear coupling between 
the signal modes and supports squeezed-light generation 
in mode {sp,ip) this way. On the other hand, values 
of squeeze variances for mode {sb^b) increase with in- 
creasing values of phase mismatch S s because of lower 
values of amplitudes of the backward-propagating pump 
mode pb- Phase relations affected by nonzero values of 
S s enable squeezed-light generation in mode (sf,Ib)' 

Increasing values of nonlinear phase mismatch 5 n i 
(S n i = 5f = Sb) result in greater values of squeeze vari- 
ances. This is demonstrated in Fig. [SJfor modes (sfj^f) 
and (sb^b)- This means that oscillations occurring 
along the waveguide and having their origin in nonzero 
values of nonlinear phase mismatch 5 n i effectively lower 
values of nonlinear coupling constants. Presence of such 
oscillations is visible in the behaviour of variances of 
quadratures ((A<j) 2 ) and ((Ap) 2 ) in mode (sf^f) (see 
Fig. |30 . The variances belonging to mode (sb^b) 
do not show oscillations (see Fig. Eb), because oscilla- 
tions are compensated in nonlinear interaction among 
the backward-propagating fields. Absence of such os- 
cillations in quantities characterizing squeezing in mode 
(sb^b) indicates that squeezing in this mode has its ori- 
gin only in nonlinear interaction among the backward- 
propagating fields. 

If we assume one of the forward-propagating signal 
and idler fields in vacuum state at the input (parametric 
downconversion) or both modes in vacuum states at the 
input (parametric amplification) we observe qualitatively 
the same behaviour of optical fields as discussed above. 

If the value of an incident amplitude of either the 
forward-propagating signal field or forward-propagating 
idler field is greater than that of the pump field (para- 
metric upconversion) squeezed light can be generated in 
compound modes (sf,«f), (s_f,?b), and (sb^b)- More- 
over, squeezed light can be reached also in compound 
modes (s f ,Pf) [(«f,Pf)] and {s b ,Pb) [(«b,Pb)] assum- 
ing A SF > A PF [Ai F > A PF ]. However, a strong negative 
influence of linear coupling affects squeezed-light gener- 




FIG. 5: Variances ((Ag) 2 } (solid line with A) and {(Ap) 2 } (*) 
and principal squeeze variance A (o) for mode (sf, If) (a) and 
(sb,ib) (b) in dependence on nonlinear phase mismatch 8 n i 
(S n i — 8f = 5b); L — 2 and values of the other parameters 
are the same as in Fig. 1. 



2.00 




FIG. 6: Principal squeeze variance A of mode (sf,«f) (solid 
line with*), (sf,«b) (A), (sb,ib) (o), (sf,Pf) (circles), and 
(sb,Pb) (squares) in dependence on linear coupling constant 



Ki (Ki 



Ks 



Ki); L = 2, A s 



10, Ai; 



An 



values of the other parameters are the same as in Fig. 1. 



ation in this regime (see Fig. HJ). We note that a mea- 
surement of squeezing in a compound mode containing 
the pump field would require a more general scheme of 
homodyne detection that would use two local oscillators, 
one of them having the carrying frequency of the pump 
field, the second one oscillating at signal- or idler-field 
frequency. 

An incident squeezed light in the forward-propagating 
pump mode does not support squeezed light genera- 
tion. This means that values of principal squeeze vari- 
ances cannot be lower than the minimum from the val- 



FIG. 7: Principal squeeze variance A of mode (sb, is) in de- 
pendence on squeeze parameter r S F and squeeze phase i? b f 
characterizing the incident field; L = 2 and values of the 
other parameters are the same as in Fig. 1. 



ues characterizing the incident light in a given mode and 
those obtained assuming an unsqueezed incident forward- 
propagating pump mode. 

If an incident squeezed light occurs in the forward- 
propagating signal mode, smaller values of squeeze vari- 
ances can be reached in modes (sb, »b) and (sb^f) com- 
pared to those reached with an unsqueezed input light. 
Values of the principle squeeze variances then depend on 
incident squeeze parameter r and incident squeeze phase 
$ [see Fig. [7| for mode (sb,ib)]- Values of principle 
squeeze variances of the other modes are greater than 
minimum from values characterizing the incident light 
and those appropriate for all unsqueezed incident beams. 



IV. SUB-POISSONIAN-LIGHT GENERATION 

Statistical properties of photoelectrons emitted inside 
a detector reflect statistical properties of integrated in- 
tensity W of the radiation impinging on the detector. 
Integrated intensity W of an optical field is defined as: 

W = iU, (27) 

where A means an electric-field-amplitude operator of 
an optical field. Detection process is characterized by 
normally-ordered moments of integrated intensity W: 



k = 2,3, 



(28) 



Statistical properties of photoelectrons are usually 
judged according to the value of Fano factor F n defined 
as: 



F n 



((An) 2 ) 
(n) 



1 



(29) 



Symbol n denotes the number of photoelectrons, An = 
n— (n), and AW — W— (W)jj. We assume that intensity 
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operator W%j of the compound mode is determined 
using the relation Wij = Wi+Wj, where Wi (Wj) denotes 
intensity operator of mode i (j). The condition F n > 1 
characterizes classical fields whereas values of F n lower 
than one can be reached only for nonclassical fields (sub- 
Poissonian light). Such fields have fluctuations in the 
number of photoelectrons suppressed below the classical 
limit (F n — 1 for a coherent state of laser radiation). 

In the following we pay attention to signal and idler 
fields that are assumed to have intensities at single- 
photon level. The planar nonlinear photonic-band-gap 
waveguide is pumped by a strong forward-propagating 
pump field. Classical amplitudes A SF , Ai F , A Sb , and 
A iB are zero in this case. Operator amplitudes A of the 
signal and idler fields are then given just by their linear 
operator corections SA. The quantity 10 Vm _1 is used 
as a unit for linear amplitude corrections. Mean values 
of intensities using this unit are then directly equal to 
mean photon numbers. 

A. Weak-interaction approximation 

The expression for Fano factor F n in Eq. (|29J) shows 
that sub-Poissonian light is generated provided that 
((AW) 2 )jv < 0. This condition cannot be fulfilled in 
a single-mode case as is indicated by the following ex- 
pressions valid in weak-interaction approximation and 
assuming coherent states with amplitudes £j for incident 
quantum linear corrections 8Aj\ 

((AW SF fU = 2\I PF \ 2 \£ SF \ 2 , 

((AW, B )V = 2|/ PB | 2 |6 B | 2 . (30) 

Constants I PF and I PB are defined in Eq. Q24[). How- 
ever, we arrive at the following expressions for compound 
modes under the same conditions: 

((A(W SF +W lF )) 2 U = 4R*{lp,C,ef,} 

+ 2|V| 2 (l + 3|6 F | 2 + 3|e iF | 2 ) 

+ 4Re {i PF i;cj: F + i PF ^: B c F } , 

((A(W SF + W SB )) 2 U = 

2(|/ PP | 2 I& F | 2 + |/ Pfl | 2 l& fl | 2 ), 

((A(W SF +W lB )) 2 U = 

2(\i PF n SF \ 2 + \i PB \ 2 \i lB \ 2 ) 

+ 4Re{-I itPF £M B + I PB ,sCM B } . 
((A(W SB +W iB )) 2 U = 4Re{/ PB C B £U 
+ 2|/ p J 2 (l + 3|6 B | 2 + 3|e i J 2 ) 

+ me{-i PB i s e SF a; B - i pb i*z* f c b } ; 

(31) 

Is = [ dzK s (z), 
Jo 

h = I dzK*(z); (32) 
Jo 



Ii, PF and Ip BlS are defined in Eq. I|26|l . Analysis of ex- 
pressions in Eq. (|31|l leads to the conclusion that negative 
variances ((AW) 2 )j^ of integrated intensity can occur in 
modes (sf^f) and (ss,is) provided that phases of the 
interacting fields are suitably chosen. The best condi- 
tions for sub-Poissonian- light generation in mode (sf, if) 
[(sb,«b)] occur provided that a,rg(I PF ^* F S,* F ) = 7r + 2ttI 
[a.Tg(I PB £t B £* B ) = 7T + 27rZ], / e N. The occurrence of 
negative variances ((AW) 2 )j^ of integrated intensity in 
mode (sf,*b) requires a stronger linear interaction than 
the nonlinear one and again a suitably chosen phases of 
the interacting optical fields as can be deduced from the 
corresponding expression in Eq. I|31|l . 



B. Numerical analysis of sub-Poissonian-light 
generation 

Complete analysis of sub-Poissonian behaviour of the 
interacting modes can be done only numerically. Pro- 
vided that states of the modes can be described in the 
framework of generalized superposition of signal and 
noise 10], normally ordered moments of integrated in- 
tensity as well as photon-number distribution can be 
determined from parameters characterizing the state in 
terms of Laguerre polynomials. Details can be found in 

Numerical analysis has shown that sub-Poissonian 
light can be generated only in compound modes (sp, If), 
(sf^b), and (ss,is). Moreover, sub-Poissonian-light 
generation requires nonzero values of the incident weak 
forward-propagating signal and idler fields. As is demon- 
strated in Fig. [SI even fields containing less than one 
photon on average are sufficient to stimulate generation 
of sub-Poissonian light. Then we can have light with val- 
ues of Fano factor F n around 0.8 in mode (sfiIf)- Scat- 
tering of signal and idler fields leads to sub-Poissonian 
statistics in modes (sf,*b) and (sb,?b). This is remark- 
able especially for mode (sb^b) being in vacuum state 
at the input (see Fig. |SJ:). The larger the incident signal- 
and idler-field amplitudes, the smaller the values of Fano 
factor F n . 

The possibility to generate sub-Poissonian light de- 
pends on length L of the waveguide (see Fig. |5J. 
Sub-Poissonian statistics of light in modes {sf,if) and 
(sfiIb) occur only for shorter waveguides. On the 
other hand, greater values of length L are convenient 
for sub-Poissonian statistics of light in mode (sb,Ib)- 
Quantities characterizing this mode composed of only 
backward-propagating fields show saturation in depen- 
dence on length L being typical for systems containing 
backward-propagating fields. 

As is indicated by analytical expressions in Eq. I|31|) 
valid in weak-interaction approximation an efficient gen- 
eration of sub-Poissonian light requires a suitable choice 
of phases of the forward-propagating fields. A strong in- 
fluence of phase ip SF of the incident forward-propagating 
signal field on values of Fano factor F n of modes (sf, «f), 
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a) 



b) 



c) 




FIG. 8: Fano factor F n for mode (sf,If) (a), (sf,ib) (b), 
and (sb,ib) (c) in dependence on amplitudes £ ap and 
linear operator amplitude corrections 8 A axe assumed to be in 
coherent states with amplitudes £ at the input; Kf = Kb = 
5 x 1(T 2 , K s = Ki = 5, L = 0.2, S s = Si = 0, S F = 8 B = 



0, ^.pj? 
arg(?s F ) 



10, A sp 
: argfe F ) 



o, €pp 







-4, 

: 0. 



0; 



(sf,ib), and (sb^b) is visualized in Fig. 

In order to reach smaller values of Fano factor F n 
greater values of the pump field A PF are necessary. The 
larger the pump amplitude A PF , the smaller the values 
of Fano factor F n . For example, assuming mode (sb^b) 

0.8 for A PF = 10 



F n ss 0.9 for A PF = 5, whereas F n 



(see Fig. [TTll- 
Linear coupling of modes influences sub-Poissonian- 
light generation as follows. The greater the values of 
linear coupling constants K s and Ki, the greater values 
of Fano factor F n are reached in mode (sF,ip). How- 
ever, even larger values of K s and Ki enable to generate 



1.1 



1.0 




^> 
1.0 1.5 

L 



FIG. 9: Fano factor F„ of mode (sf,if) (solid line with *), 
{sF,i B ) (A), and (sb, is) (o) as a function of length L; £ 3F = 
—1, £ ip = 1, and values of the other parameters are the same 
as in Fig. 8. 




FIG. 10: Fano factor F n of mode (sf,If) (solid line with *), 
(sf,ib) (A), and (sB,i B ) (o) in dependence on the phase ip SF 
[ip 3F = arg(£ SF )] in units of n; |£ SF | = 1, £ ip = 1, and values 
of the other parameters are the same as in Fig. 8. 



light with sub-Poissonian statistics. In mode (sf,«s), 
greater values of K s result in greater values of F n . On 
the other hand, greater values of Ki lead to lower values 
of F n owing to a stronger scattering between the forward- 
and backward-propagating idler modes. The behaviour 
of mode (ss,zs) is shown in Fig. ^| If strength of 
scattering between the signal and idler fields is more-less 
balanced (K s Ki), larger values of linear coupling con- 
stants K s and Ki provide smaller values of Fano factor 
F n . Otherwise greater values of K s and Ki neednot nec- 
essarily lead to smaller values of Fano factor F n . 

Complete phase matching of all interactions cannot 
be usually reached in real nonlinear photonic-band-gap 
waveguides. The effect of linear signal-field phase mis- 
match S s on values of Fano factor F n is relatively weak 
(see Fig. ^J. Greater values of 8 S decrease values of 
Fano factor F n of mode (sf^f) but they increase val- 
ues of Fano factor F n of mode (s.b,z_b). This behaviour 
can be understood claiming that nonzero values of linear 
phase mismatch 5 S effectively weaken the linear coupling 
constant between the signal modes. 

Nonzero values of nonlinear phase mismatches Sf and 
Sb effectively weaken nonlinear interaction in a waveg- 
uide and this acts against sub-Poissonian-light genera- 
tion. As is shown in Fig. ^| sub-Poissonian light cannot 
be generated for greater values of 8 n i (S 
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FIG. 14: Fano factor F n of mode (sf,If) (solid line with 
*), (sf,ib) (A), and (sb,ib) (o) in dependence on nonlinear 



phase mismatch 5 n i (S n 



5b); Us 



-1, & 



1, and 



values of the other parameters are the same as in Fig. 8. 



FIG. 11: Fano factor F n of mode (ss,is) in dependence on 
incident pump amplitude A PF and length L; £ 3F = — 1, £i p — 
1, and values of the other parameters are the same as in Fig. 
8. 





FIG. 15: Fano factor F„ of mode (sf,*f) in dependence on 
coherent amplitude £ SF and mean number of incident noisy 
photons n chiSF ; arg(£ SF ) = 0„ £ iF = —1, and values of the 
other parameters are the same as in Fig. 8. 



FIG. 12: Fano factor F„ of mode (ss,is) in dependence on 
linear coupling constants K s and Ki\ L — 0.5, £ 3F = — 1, 
£i F = 1, and values of the other parameters are the same as 
in Fig. 8. 
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FIG. 13: Fano factor F n of mode (sf,If) (solid line with 
(sf,ib) (A), and {sb,ib) (o) in dependence on linear phase 
mismatch <5 S ; f SF = — 1, ^i F = 1, and values of the other 
parameters are the same as in Fig. 8. 



all. 

Noise present in incident fields leads to suppression of 
sub-Poissonian behaviour. Generation of sub-Poissonian 
light is rather sensitive to values of incident noise, as is 
shown for mode (sf,*f) in Fig. 1151 

A nonlinear planar waveguide can be also used for the 
suppression of incident noise. This effect has its origin 
in sensitivity of the nonlinear process to the phase of an 
incident light. To be more specific, the amplification co- 
efficient of the incident light depends on its initial phase. 
If the central phase of the incident field (corresponding 
to a coherent signal amplitude) has the strongest ampli- 
fication then the noisy part (with a blurred phase) is less 
amplified on average and signal-to-noise ratio increases. 
Reduction of the incident noise is demonstrated in Fig. 
I16l for single mode sp using second reduced moment Rw 
of integrated intensity W [Rw = {W 2 )jv/ (W)jj-]. Reduc- 
tion of the incident noise occurs only for shorter lengths L 
and values of incident phases in a certain region (amplifi- 
cation of the incident field has to occur in the waveguide). 
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motion: 



FIG. 16: Second reduced moment Rw of integrated intensity 
W of mode sf in dependence on length L of the waveguide 
and incident phase ip ap [<p SF = arg(£ sp )]; \£ 3p \ = 2„ £ ip = 2, 
n c h,s F = 2, and values of the other parameters are the same 
as in Fig. 8. 



V. CONCLUSIONS 

We have analyzed optical parametric process occur- 
ring in a nonlinear planar photonic-band-gap waveg- 
uide. It has been shown that squeezed light and light 
with sub-Poissonian photon-number statistics can be 
generated only in compound modes containing a) sig- 
nal and idler forward-propagating fields, b) signal (idler) 
forward-propagating field and idler (signal) backward- 
propagating field, and c) signal and idler backward- 
propagating fields. The role of parameters characterizing 
a real nonlinear photonic-band-gap waveguide in gener- 
ation of light with nonclassical properties has been an- 
alyzed in detail and understood. Also the influence of 
incident squeezed light and incident noise has been elu- 
cidated. The waveguide can be used for the suppression 
of incident noise. The obtained analysis will be used for 
the construction of a real nonlinear photonic-band-gap 
waveguide. 
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APPENDIX A: QUANTUM DERIVATION OF 
THE NONLINEAR EQUATIONS 

Quantum description of nonlinearly interacting opti- 
cal modes is based on the construction of momentum 
operator G(z) which determines Heisenberg equations of 



dX 

dz 



G, X 



(Al) 



X stands for an arbitrary operator and [ , ] means a 
commutator. 

If a nonlinear interaction involves counterpropagating 
beams, no momentum operator can be straightforwardly 
assigned to Heisenberg equations as they are written in 
Eq. - However, we can proceed as follows QJ|. We 
assume a nonlinear interaction among all involved fields 
as if they propagate in one direction and then we can 
write the momentum operator G(z) in the form: 



a—s F ,i F ,pf 



a=s B ,ib ,Pb 



TiK s exp(iSiz)al P a SB + hK l exp{i5iz)d\ F d iB + h.c. 

2iTiK F d PF d\ F d\ F + 2iTiK B d PB d\ B d\ B + h.c. 

(A2) 

We now substitute the creation operators (al) of the 
backward-propagating fields by newly introduced fictions 
annihilation operators (b) and vice versa, i.e. 



b a) 



U PB ' 

SI 



"Sfl S B 1 ">i 

®sb *~ b\ B> Oi B *— b ig) u P] . 
Heisenberg equations in Eq. i|Alfl then have the form 
da x „ 



(A3) 



dz 

d<ii F 
dz 

d%B 
dz 

db\ 

IB 

dz 



= i{k SF ) z a SF +iK s exp(iSiz)bl B 

+ 2K F a PF a\ F , 
= i{ki P ) z d ip +iKiexp{iSiz)b} B 
+ 2K F d PF al F , 

= -i(k SB ) z b' SB - iK* exp(-iSiz)d SF 
-2K B bl B b lB , 

= -i(h B ) z bl B - iK* exp(-i6iz)d iF 
-2K R bLb 



PB S B J 



da 



IT 



dz 
dz 



IK p(L SF di F , 



^ ( F ) z @>p 
-i{k PB ) z bl B +2K B b\ B b\ 



(A4) 



Returning to the original operators a) , a in Eq. (|A4|) us- 
ing substitution in Eq. (| A3|> and transforming Eq. (|A4|) 
into interaction picture (d a (z) = A a (z) exp[z(/c a ) z z]) we 
arrive at the system of equations written in Eq. l|rj)l. 
We only have electric-field-amplitude operators instead 
of classical electric- field amplitudes occurring in Eq. JfJJl . 
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Quantum interpretation of the nonlinear process of 
three mode interaction at single-photon level then im- 
mediately provides the following conservation law of the 
overall number of photons in the interaction: 

~flZ (^s F ^s F + A\ p A iF + 2A pF A PF 
~A\ b A Sb — A\ B Ai B — 2A^ B A PB \ = 0. (A5) 

The symbol ( ) denotes mean quantum statistical value. 



APPENDIX B: IDENTITIES OBEYED BY A 
SOLUTION OF EQ. (ITT1) 

Solution of Eq. Ullfl for the annihilation operators of 
linear operator corrections can be expressed as: 



( SA(L) 
\6BHL) 



«n «i2 \ ( 5A(0) \ 
u 2 i u 22 ) \6Bi(0) J 



V21 V 2 2 



and 



5A{z) 



'5A SF {z)' 
SA iF (z) 



SB\z) 



/<5it(0) 
V SB(0) 



'5A SB {z)' 



(Bl) 



(B2) 



Quantum method of the derivation of Eq. © de- 
scribed in Appendix A indicates that the following com- 
mutation relations have to be fulfilled: 



[5A t (L),SA k (L)} = 0, 
[SML),SAl(L)] = 6 it 
[5Ai(L),SB k (L)} = 0, 



[SAi(L), 5B\(L)] 
[6Bj(L),SB k (L)] 
[5b\{L),5B{{L)} 



0. 
0. 



(B3) 



Substitution of expressions in Eq. IjBljl into relations 
in Eq. I|B3(1 provides identities for the matrices u and v 
defined in Eq. (|B1|I : 

[(un)ij(vu)kj - (ui2)ij(vi2)kj ~ (vu)ij(un)kj 

3 

+ (v 1 2)tj(ui 2 )kj] = 0, 

3 

+ ( u 12)y(«* 2 )fej] = $ik, 

Y [( M ll)^"( U 2l)fci - (Ul2)ij(u2 2 )kj - («ll)ii(«2l)fcj 
3 

+ ( u 12)y(«22)fcj] = 0, 

Y [(Un)ij(v2l)kj ~ {Ul2)ij(V22)kj ~ (Vll)ij(U2l) k j 
3 

+ {vi2)ij(u22)kj] = 0, 

Y [( M 2l) 4 i(u2l)fci - (U22)ij(l4 2 )kj - («2l)ii(«2l)fcj 
3 

+ (v 2 2)ij{V2* 2 )kj] = -$ik, 

^ [(U2l)ij(v 2 l)kj - (u22)ij(V22)kj ~ {V2l)ij(U2l)kj 
3 

+ (V22)ij(u 2 2)kj] = 0. 

(B4) 

Identities given in Eq. (|B4|) are useful for controlling 
precision of numerical integration of the corresponding 
differential equations. 
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